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29 ■ Abstract 
0\ ■ 

, Direct dark matter detection is considered in the Next-to-Minimal Su- 

\ per symmetric Standard Model (NMSSM). The effective neutrahno-quark La- 

grangian is obtained and event rates are calculated for the ''^Ge isotope. Ac- 
^0 ■ celerator and cosmological constraints on the NMSSM parameter space are 

\ included. By means of scanning the parameter space at the Fermi scale we 

show that the lightest neutralino could be detected in dark matter experiments 
. with sizable event rate. 

; PACS number(s): 14.80.Ly, 95.35+d, 98.80Cq 
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0\ ; 1 Introduction 

^ ; 

Ph, In not too far future new very sensitive dark matter (DM) detectors [|T1^0 may 

start to operate, and one expects new, very important data from these experiments. 
^ \ The future experimental progress forces investigators to know better the variety and 

property of the dark matter particles. The lightest supersymmetric particle (LSP), 



X 



the neutralino, is considered now as a most promising candidate, which may compose 
the main fraction of so-called cold dark matter. The prospects of the direct and 
5^ \ indirect detections of the LSP have comprehensively been investigated in the 

various versions of the Minimal Supersymmetric Standard Model (MSSM) 0. 

In this paper we consider direct detection of this relic LSP in the Next-To-Minimal 
Supersymmetric Standard Model (NMSSM) @, 0]. The Higgs sector of the NMSSM 
contains five physical neutral Higgs bosons, three Higgs scalars, two pseudoscalars, 
and two degenerate physical charged Higgs particles C^. The neutralino sector is 
extended to five neutralinos instead of four in the MSSM. The remaining particle 
content is identical with that of the MSSM. 

The NMSSM is mainly motivated by its potential to eliminate the so-called /i 
problem of the MSSM p], where the origin of the parameter in the superpotential 
W^MSSM = is not understood. For phenomenological reasons it has to be of 
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the order of the electroweak scale, while the "natural" mass scale would be of the 
order of the GUT or Planck scale. This problem is evaded in the NMSSM where 
the n term in the superpotential is dynamically generated through fi = Xx with a 
dimensionless coupling A and the vacuum expectation value x of the Higgs singlet. 
Another essential feature of the NMSSM is the fact that the mass bounds for the 
Higgs bosons and neutralinos are weakened. While in the MSSM experimental data 
imply a lower mass bound of about 20 GeV for the LSP f^, very light or massless 



neutralinos and Higgs bosons are not excluded in the NMSSM |10, 11 1. Furthermore 



the upper tree level mass bound for the lightest Higgs scalar of the MSSM 

ml < m| cos^ 2(3 (1) 

is increased to < m| cos^ 2/9 + X^^vf + f 2) sin^ 2/?. Taking into account the weak 
coupling of the singlet Higgs the NMSSM may still remain a viable model when the 
MSSM can be ruled out due to (|l]). 

The above arguments make an intensive study of the NMSSM phenomenology very 
desirable. Previously the Higgs and neutralino sectors of the NMSSM were carefully 
studied in |TD[-|[T3|. The calculation of the LSP relic abundance in the NMSSM was 



performed for the first time in |]T6| and recently in [|T] 

The outline of this paper is as follows. In Sec. ^ we describe the Lagrangian of 
the NMSSM. Since the additional singlet superfield of the NMSSM leads to extended 
Higgs and neutralino sectors, we present the Higgs and neutralino mixings. Section |^ 
collects formulas relevant for calculation of the event rate for direct dark matter 
detection in the framework of the NMSSM. In Sec. ^ we discuss the constraints on 
the NMSSM parameter space which are used in our analysis. In Sec. ^ we shortly 
describe our numerical procedure and discuss the results obtained. Sec. |] contains a 
conclusion. 



2 The Lagrangian of the NMSSM 

The NMSSM superpotential is [Q {eu = -621 = 1): 

W = XeijHiHiN - -kN^ 

3 

+KeijQVHi - h^e^.QWHi - Ke.fVRHl (2) 

where Hi = {Hf,H^) and H2 = {H^^H^) are the SU(2) Higgs doublets with hyper- 
charge —1/2 and 1/2 and is the Higgs singlet with hypercharge 0. The notation 
of the fermion doublets and singlets is conventional, generation indices are omitted. 
Contrary to the MSSM, the superpotential of the NMSSM consists only of trilinear 
terms with dimensionless couplings. 

The electroweak gauge-symmetry SU{2)i x f/(l)y is spontaneously broken to the 
electromagnetic gauge-symmetry U{l)em by the Higgs VEVs = Vi with i = 1,2 



and (A^) = x, where v = Jvl + f f = 174 GeV, tan/? = V2/V1. 
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The most general supersymmetry breaking potential can be written as 

V^oft = ml\H,\' + ml\H,\' + ml\N\' 

+ml\Q\^ + ml\U\^ + ml\b\^ + ml\L\^ + ml\R\^ 

-{XAxe,,HlHiN + h.c.) - {hA,N^ + h.c.) 

+iKAueijQ'UHi - hdADSi.Q'DHi - KAESifURHi + h.c.) 

+^MA"A" + ^M'A'A'. (3) 

As free parameters appear the ratio of the doublet vacuum expectation values, 
tan (3, the singlet vacuum expectation value x, the couplings in the superpotential 
A and fc, the parameters A\, Ak, as well as Au, Ad, Ae (for three generations) in 
the supersymmetry breaking potential, the gaugino mass parameters M and M', and 
the scalar mass parameters for the Higgs bosons mi^2,3, squarks ruQ^u^D and sleptons 

The minimization conditions for the scalar potential dV/dvi^2 = 0, dV/dx = 
eliminate three parameters of the Higgs sector which are normally chosen to be ml, 
m\ and m| . Then at tree level the elements of the symmetric CP-even mass squared 
matrix = (M-"^^) become in the basis {Hi, H2, N) 

Mfi^ = ^vl{g'^ + g^) + Xxtanp{Ax + kx), 
Mf/ = -Xx{Ax + kx)+ViV2{2X'-^g''-^g') 



• 2 

'13 

rS'2 _ 1„,2/ /2 , 2x 



= 2X'^ViX - 2Xkxv2 - XAxV2, 



^22 = -viig" + g') + Xxcot(3{Ax + kx), 
Mfg^ = 2X'^V2X -2Xkxvi - XAxVi, 
= Ak'^x^ - kAkX + 



X 

2 



In the same way one finds for the elements of the CP-odd matrix M.p 



M[i = Xx{Ax + kx) tan f3, M^^ = Xx{Ax + kx), 
Mi^^ = Xv2{Ax - 2kx), M^2 = >^x{Ax + kx) cot 
M|;^ = Xvi{Ax - 2kx), M^^ = XAx^ + 4:XkviV2 + 3kAkX, 



X 

and for the charged Higgs matrix one obtains 
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In our numerical analysis we have included 1-loop radiative corrections to Higgs mass 
matrices following [|14[ |l^. 

Assuming CP conservation in the Higgs sector, the Higgs matrices are diagonalized 
by the real orthogonal 3x3 matrices and U^, respectively, 



Diag(m^^,m^^,0) = U^^M^U 



where < < "^Ss and mp^ < mp^ denote the masses of the mass eigenstates 
of the neutral scalar Higgs bosons Sa (a = 1,2,3) and neutral pseudoscalar Higgs 



bosons Pa (a = 1,2) p 



With fixed parameters of the Higgs sector the masses and mixings of the neutrali- 
nos are determined by the two further parameters M and M' of the Lagrangian 

In this basis the symmetric mass matrix M of the neutralinos has the form: 

/ M' —mz sinOyy cos/3 mz sinOw sin/? \ 

M mz cos^iy cos/3 —mz cos 9w sin/3 

—mz smOw cos/3 mz cos^iy cos/3 Ax A 1^2 

mz sin 0(4/ sin/3 —mz cos 9w sin/3 Ax A t;i 

V \V2 Xvi -2k X J 

The mass of the neutralinos is obtained by diagonalizing the mass matrix M with 
the orthogonal matrix A^: 

with x° = ATij^j and Mjiag = N M N^. 

The neutralinos {i = 1-5) are ordered with increasing mass thus x = Xi 
the LSP neutralino. The matrix elements Afij = 1-5) describe the composition 
of the neutralino Xi the basis For example the bino fraction of the lightest 
neutralino is given by A/'fi and the singlino fraction of this neutralino by A/^s- 



3 Neutralino-nucleus elastic scattering 

A dark matter event is elastic scattering of a DM neutralino from a target nucleus 
producing a nuclear recoil which can be detected by a suitable detector. The corre- 
sponding event rate depends on the distribution of the DM neutralinos in the solar 
vicinity and the cross section of neutralino-nucleus elastic scattering. 
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The relevant low-energy effective neutralino-quark Lagrangian can be written in 

(4) 



a general form as |J, |1^, [2^, pi 

^eff 



m„ 



w 



m 



4 / ' 



where terms with vector and pseudoscalar quark currents are omitted being negligible 
in the case of non-relativistic DM neutralinos with typical velocities v ~ lO^^c. 
The coefficients in the effective Lagrangian (H) have the form: 



Here 









L 2 






m?^ - 









ml ( 1 



cos^ 05l + sin^ Oq (pl^) 
sin^ Oq (pl^ + cos^ ^5 (plji) 



ip2 



+ 



9l 



(m^ + mg)2 m?2 - (m^ + 



Mh/ Pq sin20g T3 (A/Is - tan ^H/A/"! 



11; 



92 



_ d 
4 

+ n 



1 



a=l,2,3 

'COS^ 9q (pqL - Sin^ 9q (pqR COS^ 6'g </)g/J - Sin^ 9q (pqL 



m 



+ sm2eq{ 



gi 
rrir, 



1 



Mw , 



m 



X 



m?i - (m^ + mqY - (m^ + ruq) 



V, 



aq 



,1 1 



•^al 



cos/5 



Qan = (M2-tanVA/'ii)[t/i'iA/'i3-f/i'2M 



14 



= Afun+AfniQ-n)tanew, 
= tanOw Q A/11, 



15; 



-1 ^xA/'l4 ,1 

■2+^^W + ^2 



cos P 



(5) 



(6) 
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The coefficients Aq and Cq take into account squark mixing — q^ and the 
contributions of all CP-even Higgs bosons. Under the assumption A = A; = these 



formulas coincide with the relevant formulas in the MSSM [IS]. 

A general representation of the differential cross section of neutralino-nucleus scat- 
tering can be given in terms of three spin-dependent Tij{q^) and one spin-independent 



Ts{q ) form factors as follows p2| 



^{v, q^) = ^(flo ■ -^00 (?^) + «oai ■ J^io(g^) 

+ al.J^W) + cl.A^J^Iiq^)). (7) 



The last term corresponding to the spin-independent scalar interaction gains coher- 
ent enhancement A"^ {A is the atomic weight of the nucleus in the reaction). The 
coefficients ao,i,Co do not depend on nuclear structure and relate to the parameters 
Aq, Cq of the effective Lagrangian (^) and to parameters characterizing the nucleon 



structure. In what follows we use notations and definitions of our paper [23|. 

An experimentally observable quantity is the differential event rate per unit mass 
of the target material 



dR 

IK. 





/ 




Jv 



dvf{v)v%{v,Er), q^ = 2MAE, 



dq 



Here f{v) is the velocity distribution of neutralinos in the earth's frame which is 
usually assumed to be a Maxwellian distribution in the galactic frame. is the 
number density of the target nuclei. Vmax = Vesc ~ 600 km/s and = 0.3 GeV-cm~^ 
are the escape velocity and the mass density of the relic neutralinos in the solar 

1 /2 

vicinity; Vmin = {MAEr/2Ml^^^) with Ma and Mred being the mass of nucleus A 
and the reduced mass of the neutralino-nucleus system, respectively. 

The differential event rate is the most appropriate quantity for comparing with 
the observed recoil spectrum and allows one to take properly into account spectral 
characteristics of a specific detector and to separate the background. However, in 
many cases the total event rate R integrated over the wholere kinematic domain of 
the recoil energy is sufficient. It is widely employed in theoretical papers for estimating 
the prospects for DM detection, ignoring experimental complications which may occur 
on the way. In the present paper we are going to perform a general analysis aimed 
at searching for domains with large values of the event rate R like those reported in 



24 . This is the reason why we use in the analysis the total event rate R. 



4 Constraints on the NMSSM parameter space 

Assuming that the neutralinos form a dominant part of the DM in the universe one 
obtains a cosmological constraint on the neutralino relic density. The present lifetime 
of the universe is at least lO^'' years, which implies an upper limit on the expansion 
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rate and correspondingly on the total relic abundance. Assuming ho > OA one finds 
that the contribution of each relic particle species x has to obey pSf: fix^o < 1; 



where the relic density parameter il^. = Px/ Pc is the ratio of the relic neutralino mass 
density to the critical one pc = 1.88 ■ lO'^^hpg-cm"^. 

We calculate fi^^/ig following the standard procedure on the basis of the approxi- 



mate formula [26, 27 



QX = 2.13xl0-"f^l (JI:l_)\n'/( ^'V! , V (8) 
° \T^J \2.7K"J ^ \axF + bxl/2j ^ ' 

Here is the present day photon temperature, T^/T^ is the reheating factor, xp = 
Tp/rrix ~ 1/20, Tp is the neutralino freeze-out temperature, and Np is the total 
number of degrees of freedom at Tp. The coefficients a, h are determined from the 
non-relativistic expansion (aann.v) ~ a + 6x of the thermally averaged cross section 
of neutralino annihilation in the NMSSM. We adopt an approximate treatment not 
taking into account complications, which occur when the expansion fails [|28[. We 
take into account all possible channels of the X'X annihilation. The complete list of 
the relevant formulas in the NMSSM can be found in ITT . 



Since the neutralinos are mixtures of gauginos, higgsinos, and singlino the an- 
nihilation can occur both, via s-channel exchange of the Z° and Higgs bosons and 
t-channel exchange of a scalar particle, like a selectron. This constrains the parameter 
space, as discussed by many groups pO |. 



In the analysis we ignore possible rescaling of the local neutralino density p which 
may occur in the region of the NMSSM parameter space where ^x^o < 0.025 |^ 



32| , p^ . If the neutralino is accepted as a dominant part of the DM its density has 
to exceed the quoted limiting value 0.025. Otherwise the presence of additional DM 
components should be taken into account, for instance, by the mentioned rescaling 
ansatz. However, the halo density is known to be very uncertain. Therefore, one 
can expect that the rescaling takes place in a small domain of the parameter space. 
Another point is that the SUSY solution of the DM problem with such low neutralino 
density becomes questionable. We assume neutralinos to be a dominant component 
of the DM halo of our galaxy with a density px = 0.3 GeV-cm~^ in the solar vicinity 
and disregard in the analysis points with ^x^o ^ 0.025. 

The parameter space of the NMSSM and the masses of the supersymmetric parti- 
cles are constrained by the results from the high energy colliders LEP at CERN and 
Tevatron at Fermilab |Ty, |11| . A key role for the production of Higgs bosons at e~^e~ 
colliders plays the Higgs coupling to Z bosons, while neutralino production at LEP 
crucially depends on the Zx^x^ coupling which is formally identical in NMSSM and 
MSSM and differs only by the neutralino mixing. All those couplings are suppressed 
in the NMSSM if the respective neutralinos or Higgs bosons have significant singlet 
components. Therefore NMSSM neutralino and Higgs mass bounds are much weaker 
than in the minimal model \T2]. The consequences from the negative neutralino search 



at LEP for the parameter space and the neutralino masses have been studied in [jTO|. 
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In it is shown that a very hght NMSSM neutrahno cannot even be ruled out at 



LEP2. 

We used the following constraints from LEP. For new physics contributing to the 
total Z width Ar(Z + 2' < 23 MeV. For new physics contributing 

to the invisible Z width Ar(Z Xi'Xj') < 8 MeV. From the direct neutralino search 
B{Z < 2 X 10-5 for j = 2,..., 5, and B{Z < 5 x 10-^ for 

i, j = 2, . . . ,5. The results of LEP searches for S^Z and S^Z* productions which 
impose restrictions on the SJZZ couplings were included in our analysis. We have 
also included the experimental bounds from the direct search for pseudoscalar Higgs 
bosons produced together with a Higgs scalar at LEP P, but this in accordance with 



does not significantly affect the excluded parameter domain. 
In our numerical analysis we use the following experimental restrictions for the 
SUSY particle spectrum in the NMSSM: m-+ > 90 GeV for charginos, mj> > 80 GeV 
for sneutrinos, mg^ > 80 GeV for selectrons, > 150 GeV for squarks, m^^ > 
60 GeV for light stop, raii+ > 65 GeV for charged Higgses and > 1 GeV for the 
light scalar neutral Higgs. In fact, it appeared that all above-mentioned constraints 
do not allow ms^ to be smaller then 20 GeV. 



5 Numerical Analysis 

Randomly scanned parameters of the NMSSM at the Fermi scale are the following: 
the gaugino mass parameters M' and M, the ratio of the doublet vacuum expec- 
tation values, tan/5, the singlet vacuum expectation value x, the couplings in the 
superpotential A and k, squared squark mass parameters m^^ ^ for the first two gen- 
erations and rriQ^ for the third one, the parameters A\, Ak, as well as At for the third 
generation. The parameters are varied in the intervals given below 



-1000 GeV 


< 


M' 


< 


1000 GeV 


-2000 GeV 


< 


M 


< 


2000 GeV 


1 


< 


tan/3 


< 


50 


GeV 


< 


X 


< 


10000 GeV 


-0.87 


< 


X 


< 


0.87 


-0.63 


< 


k 


< 


0.63 


100 GeV^ 


< 




< 


1000000 GeV^ 


100 GeV^ 


< 




< 


1000000 GeV^ 


-2000 GeV 


< 


At 


< 


2000 GeV 


-2000 GeV 


< 


Ax 


< 


2000 GeV 


-2000 GeV 


< 


Ak 


< 


2000 GeV. 



For simplicity the other sfermion mass parameters fnfj_^^, Tn'j-)_^^, m|;^ and m|,, 

m^, "H^ig, ""^Bg, are chosen to be equal to m^^ and toq^, respectively. Therefore masses 
of the sfermions in the same generation differ only due to the D-term contribution. 
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Other parameters (except At) of the supersymmetry breaking potential Au, Ad, Ae 
(for all three generations) are fixed to be zero. 

The main results of our scan are presented in Fig. 1 in the form of scatter plots. 
Given in Fig. 1 are the total event rate R for '^^Ge, and the LSP gaugino fraction 
{■Mil +-^1^2)) singlino fraction (A/'is), and finally relic density parameter fl^hl versus 
the LSP mass. The left panel in Fig. 1 presents the above-mentioned observables 
obtained without taking into account the cosmological relic density constraint. 

In this case the total expected event rate R reaches values up to about 50 events 
per day and per 1 kg of the ^"^Ge isotope. As on can see from Fig. 1 the small-mass 
LSP (less then about 100 GeV) are mostly gauginos, with very small admixture of 
the singlino component. Large masses of the LSP (larger then 100 GeV) correspond 
to sizable gaugino and singlino components together perhaps with some higgsino 
fraction. 

The results of implementation of the cosmological constraint 

0.025 < n^hl < 1. 

one can see in the right panel of Fig. 1. There is approximately a 5- fold reduction of 
the number of the points which fulfill all restrictions in this case. Nevertheless quite 
large values of event rate R (above 1 event/day/kg) still survive the cosmological 
constraint. The lower bound for the mass of the LSP now becomes about 3-5 GeV. 
The gaugino component becomes more significant, but the singlino fraction can not be 
completely ruled out especially for large masses of the LSP. The higgsino component 
of the LSP remains still possible only for LSP masses in the vicinity of Mz- 

For illustration in Fig. 2 we present the calculated event rate R as function of the 
mass of the lightest scalar Higgs boson, rug^. The largest values of R are concentrated 
mostly in the region of quite large masses rrig^, where LEP constraints are not very 
significant. The upper bound for m^^ is also clear seen. 

6 Conclusion 

In the paper we address the question whether the Next-to-Minimal Supersymmetric 
Standard Model can be attractive from the point of view of the direct detection of 
neutralinos provided the neutralino is the stable LSP. 

To answer the question we derived the effective low energy neutralino-quark La- 
grangian, which takes into account the contributions of extra scalar Higgs boson and 
extra neutralino. On this basis we calculated the total direct-dark- matter-detection 
event rate in ^^Ge as a representative isotope which is interesting for construction of 
a realistic dark matter detector. We analyzed the NMSSM taking into account the 
available accelerator and cosmological constraints by means of random scan of the 
NMSSM parameter space at the Fermi scale. We demonstrated that the cosmological 
constraint does not rule out domains in the parameter space which correspond to 
quite sizable event rate in a germanium detector. 
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Due to relaxation of the gaugino unification condition, contrary to previous con- 



sideration we found domains in the parameter space where hghtest neutrahnos 
have quite small masses (about 3 GeV), acceptable relic abundance and sufficiently 
large expected event rate for direct detection with a ''^Ge-detector. 

Therefore the NMSSM looks not worse then the MSSM from the point of view of 
direct dark matter detection. The question arises: Is it possible to distinguish MSSM 
and NMSSM by means of direct dark matter detection of LSP? It is a problem to 
be solved in future. The question can disappear by itself if negative search for light 
Higgs with LHC rules out the MSSM. As already mentioned in the introduction 
the NMSSM can bypass the most crucial constraint for the MSSM with the upper 
bound for the light Higgs boson ([l|). Therefore the NMSSM might remain a viable 
theoretical background for direct dark matter search for relic neutralinos in the post- 
MSSM epoch. 

We thank S.G. Kovalenko for helpful discussions. The investigation was sup- 
ported in part (V.A.B.) by Grant GNTP 215 from Russian Ministry of Science and 
by joint Grant 96-02-00082 from Russian Foundation for Basic Research and Deutsche 
Forschungsgemeinschaft . 
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Figure captions 

Fig. 1: Total event rate R for ^^Ge, the LSP gaugino fraction (A/'fi + ■^12)) singlino 
fraction (jVfs), and the relic abundance parameter f^^^o versus the LSP mass (from 
up to down). The left (right) panel presents results obtained without (with) taking 
into account the cosmological relic density constraint. 

Fig. 2: Total event rate R for ''^Ge as function of the mass of the lightest scalar 
Higgs boson rrig^. 
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